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ABSTRACT 



SuperWASP light curves for 53 WUMa-type eclipsing binary (EB) candidates, identified in previous work as being close to the con- 
tact binary short-period limit, were studied for evidence of period change. The orbital periods of most of the stars were confirmed, and 
period decrease, significant at more than 5 cr, was observed in three objects: 1SWASP J1743 10.98+432709.6 (-0.055 ± 0.003 s yr" 1 ), 
1SWASP J133105.91 + 121538.0 (-0.075 ±0.013 s yr 1 ) and 1SWASP J234401. 81-212229.1 (-0.313 + 0.019 s yr 1 ). The magnitudes 
of the observed period changes cannot be explained by magnetic braking or gravitational radiation effects, and are most likely pri- 
marily due to unstable mass transfer from primary to secondary components, possibly accompanied by unstable mass and angular 
momentum loss from the systems. If these period decreases persist, the systems could merge on a relatively short timescale. 
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1. Introduction 

The SuperWASP project has been surveying bright (V ~8- 
15 mag) stars across almost the whole sky since 2004, look- 
ing f or photometric vari ations indicative of exoplanetary tran- 
sits (Pollacco et al. 2006). However, in addition to 67 exoplanets 
announced by early February 2012Q, it has discovered several 
tens of thousands of new variable stars, the majority of which 
appear to be eclipsing binaries (EB) (Payne in prep.). The cam- 
eras' observational cadence permits detection of periods from 
several months do wn to about 20 minu tes (e.g. pulsations in Am 
stars measured bv lSmallev et al.ll2.Ql ll) . Hence the SuperWASP 
archive represents a useful source of evidence for short-period 
EBs. 

W UMa-type variables, which usually represent contact bi- 
naries, exhibit a fairly sharp lo wer limit to the ir observed or- 
bital periods, of around 0.22 d dRucin ski 2007). Suggested ex- 
planations for th is cut-off have involved the stars reaching full 
convective limit (Rucinski 1992) or an effective low mass limit, 
given the finite age of the Universe, on the assumption that an- 
gular momentum loss rates are related to stellar masses (IStepienl 
2006). How ever, the convec tive limit would not provide a full 
explanation (Rucinski 1992), and Stepien's low mass limit does 
not appear to ho ld true for some recently-disc overed sys tems 
(iJiang et al.ll201ll) e.g. GSC 23 14-0530. Hence lJiang et alJ have 
suggested that dynamic mass transfer instability develops in de- 
tached binary systems possessing both low initial primary mass 
and a low mass ratio, once the primary fills its Roche lobe, lead- 
ing to a rapid merger of the binary. Stable mass transfer would 
however lead to formation of a contact binary at the point of 



http://exoplanet.eu/catalog.php 



Roche lobe overflow, explaining the short-period limit for con- 
tact b inaries wi th appropriate primary mass and mass ratio. 

iNorton et al.l (1201 ll) presented evidence for 53 candidate 
WUMa-type EBs observed using SuperWASP, with periods 
close to the short-period limit, of which 48 were new discov- 
eries. Some of these, when folded at their determined period, 
failed to give well-defined light curves, which could be the re- 
sult of period variability. Therefore our aim here was to deter- 
mine whether any of the 53 candidate EBs showed evidence 
of period changes during their observation by SuperWASP. 
Period decreases might indicate a detached or semi-detached bi- 
nary undergoing unstable mass transfer likely to end in merger, 
or even a contact b inary approaching merger (as observed by 
iTvlenda et alJl201lh . 



2. Method 

The best orbital period values for the 53 stars were checked 
using a three-step custom-written IDL procedure. First, an ap- 
proximation to the half -period was found by fitting a sinusoidal 
funct i on, using the Lev enberg-Marquardt algorithm dLevenbergl 
119441: iMarauardil 19631) . to light curve data from three optimal 
nights of observation. The inverse-variance weighted average 
half -period, and its uncertainties, were then used to constrain 
the search range for a Lomb-Scargle p eriodogram dLomblll976t 
IScargld 19821 Home & Baliunasfl 986) to determine a more pre- 
cise estimate for the half-period, which was doubled to give an 
estimate for the orbital period. Finally, the light curve data set 
was folded on trial periods within a narrow range, and the dis- 
persion of flux values within phase bins was minimized to find 
the optimum period to the nearest 0.01 s. 
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A reference light curve minimum was then automatically se- 
lected as near to the middle of the data set as possible (since the 
period here would be expected to approximate the 'average' pe- 
riod found for the whole data set in cases where the period was 
changing linearly). The time of zero-point minimum to was more 
precisely found by fitting a sinusoidal function to nearby points; 
the optimum period length P ca j c was then used to calculate ex- 
pected times of minimum light (C = to + P CSL i c E, where E is the 
period number, or epoch, relative to to) within the full time-range 
observed. Actual times of minimum O near these values were 
determined using quadratic fitting (which proved more robust 
than Gaussian or sinusoidal fits), with poor fits and inadequately- 
sampled minima being automatically rejected. 

In this way a series of observed minus calculated (O-C) val- 
ues were obtained, and plotted against epoch. A second series 
was found for the other set of light curve minima (since EBs 
exhibit a primary and secondary eclipse during each orbital pe- 
riod). Such O-C diagrams will be expected to follow a quadratic 
curve if the period is changing linearly, allowing the determina- 
tion of the rate of change of period with respect to epoch. Since 
O = t + P(E)E when P = P(E), we have: 



O - C = P(E)E - P calc £ = aE 2 + bE + c, 
so, differentiating with respect to E: 



dP 
dE 



E + (P(E) - P calc ) = 2aE + b, 



(1) 



(2) 



giving dP/dE = 2a. Therefore an attempt was made to fit a 
quadratic function to each O-C series, from which a value for 
the rate of period change (per period) could be determined. The 
inverse-variance weighted average period change (across both 
O-C series) was thus found for each star, and could be converted 
into an approximate rate of change dP/dt in s yr" 1 . 

The validity of the method and code described above was 
tested using several synthetic data sets. Even with noisy data, it 
proved able to recover input periods correct to 6 to 7 significant 
figures (s.f.), and input linear period changes correct to 0. 1 s yr" 1 , 
significant at over 7 <x. The latter error is largely due to the ap- 
proximation of a continuous function of time by a discontinu- 
ous function o f epoch in the custom ary differential equation as 
pointed out bv lKopal & Kurthl dl957l) . 

3. Results 

The optimum orbital periods found for the 53 stars are given in 
TableQ] In the majo rity of cases, identic al periods (to 4 or 5 s.f.) 
to those found in Norton et al] (1201 ll) were obt ained. Object 
J0554's period was given incorrectly in table 1 of [Norton et all 
but correctly in the appendix. For nine objects, a period ap- 
proximately 1 /8 longer was obtained by this method; manual 
checks indicated that the longer periods found here were more 
plausible (e.g. the folded light curve was more clearly defined; 
the primary and secondary minima exhibited more difference 
in depths). In the earlier paper, only candidate periods below 
20000 s (~ 0.2315 d) had been investigated, allowing these 
longer periods to be missed. Therefore the 53 stars probably do 
not all possess periods below 20000 s, but since all are below 
about 22600 s (~ 0.2616 d), they still constitute a useful sample 
of notably short-period EB candidates. 

Table Q] also gives the period change results. Object J1738 
only had enough minima data for a single O-C series, so was 
excluded. Object J0930 appeared to show significant period in- 
crease, but its folded light curve and O-C diagram suggested an 
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Fig. 1. Light curve for object J 1743, folded at period of 
22300.517 s, with mean values for binned data overplotted. 
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Fig. 2. O-C diagram for object J1743, with best fit quadratic 
curves overplotted. Black stars and solid line indicate primary 
eclipse series; grey open diamonds and dashed line indicate sec- 
ondary eclipse series. Period decrease significant at 21 <x is indi- 
cated. 



anomaly in the data, and a check of the sky within 1' of the ob- 
ject revealed a nearby bright star which appears to have contam- 
inated the light curve; this object was also therefore excluded 
at this stage. Of the remaining 51 objects, 17 showed period 
change inconsistent with zero, at 1 <x confidence, with 6 cases 
of apparent period increase and 1 1 of decrease. However, given 
the number of objects studied, we might expect approximately 
this number of marginally-significant results purely by chance, 
and so we do not claim that all these 17 results indicate genuine 
period changes (although this is a possibility). 
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Table 1. Period results for 53 SuperWASP eclipsing binaries. 



SuperWASP ID 


Original 


Short 


P 


SP 


P 


SP 


P 
P 


Significance 


(Jhhmmss.ss±ddmmss.s) 


Number" 


Name 


(s) 


(s) 


(s yr 1 ) 


(s yr 1 ) 


(yr 1 ) 


(cr) 


J000437. 82+033301.2 


3 


J0004 


22593.528" 


0.009 


+0.095 


0.087 


+4.2 x 10~ b 


1 


J003033.05+574347.6 


27 


J0030 


19580.280 


0.004 










J004050.63+071613.9 


15 


J0040 


19809.204 


0.004 










J022050.85+332047.6 r 


53 


J0220 


16643.642 


0.009 










J022727.03+1 15641.7 


47 


J0227 


18226.361 


0.009 










J030749.87-365201.7 


26 


J0307 


19584.395 


0.009 










J034439.97+030425.5 


11 


J0344 


19861.395 


0.009 










J040615.79-425002.3 


34 


J0406 


19209.980 


0.009 


-0.029 


0.026 


-1.5 x 10~ 6 


1 


J041 120.40-230232.3 


43 


J0411 


18690.406 


0.009 










J041655.13-492709.8 


5 


J0416 


19959.767 


0.009 










J044132.96+440613.7 


19 


J0441 


19712.787 


0.006 










J050904.45-074144.4 


13 


J0509 


19835.305 


0.009 










J051501. 18-021948.7 


6 


J0515 


19950.121 


0.009 










J052036.84+030402.1 


1 


J0520 


19993.301 


0.009 


-0.037 


0.025 


-1.8 x 10~ 6 


1 


J055418.43+442549.8 


39 


J0554 


18878.004 


0.003 










J064501.21+342154.9 


37 


J0645 


21480.415* 


0.009 










J074658.62+224448.5 


38 


J0746 


19081.403 


0.005 










J080150.03+471433.8 


42 


J0801 


18793.167 


0.007 










J092328.76+435044.8 


50 


J0923 


20292.982* 


0.004 


+ 1.90 


0.95 


+9.37 x 10~ 5 


2 


J092756.25-391 119.2 


30 


J0927 


19469.696 


0.009 


-0.10 


0.05 


-5.1 x 10~ 6 


1 


J093012.84+533859.6 


22 


J0930 


19674.550 


0.001 


(+1.40) 


(0.13) 


(+7.11 x 10" 5 ) 




J111931.48-395048.2 f 


14 


J1119 


19827.651 


0.007 


+0.033 


0.014 


+ 1.7 x 10~ 6 


2 


J 114929.22-423049.0 


23 


J1149 


19639.554 


0.009 










Jl 15557.80+072010.8 


25 


Jl 155 


19614.222 


0.009 










Jl 15605.88-091300.5 


48 


J1156 


18222.626 


0.005 










J1201 10.98-220210.8 


24 


J1201 


19627.832 


0.009 


-0.21 


0.09 


-1.1 x 10~ 5 


2 


J121906.35-240056.9 


29 


J1219 


19558.186 


0.009 










J130920.49-340919.9 


33 


J 1309 


19253.528 


0.009 










J133105.91 + 121538.0 


41 


J1331 


18836.380 


0.004 


-0.075 


0.013 


-4.0 x 10~ 6 


5 


J142312.63-222425.1 


51 


J1423 


18112.792 


0.009 










J150822.80-054236.9 


10 


J1508 


22469.214* 


0.007 


-0.026 


0.016 


-1.2 x 10~ 6 


1 


J15 1652.90+004835.8 


49 


J1516 


18207.370 


0.009 










J155822.10-025604.8 


9 


J1558 


22470.705* 


0.009 


-0.19 


0.13 


-8.3 x 10~ 6 


1 


J160156.04+202821.6 


28 


J1601 


19572.134 


0.009 


+0.12 


0.05 


+6.0 x 10~ 6 


2 


J161334.28-284706.7 


12 


J1613 


19852.817 


0.004 


-0.063 


0.055 


-3.2 x 10~ 6 


1 


J170240.07+151123.5 / 


4 


J 1702 


22589.758* 


0.004 










J173003.21+344509.4 


32 


J1730 


19328.922 


0.009 










J173828.46+1 11 150.2 


35 


J1738 


21546.731* 


0.009 


s 








J174310.98+432709.6'' 


17 


J 1743 


22300.517* 


0.009 


-0.0546 


0.0025 


-2.45 x 10~ 6 


21 


J180947.64+490255.0' 


20 


J 1809 


19688.490 


0.009 


+0.028 


0.010 


+ 1.4x 10~ 6 


2 


J183738. 17+402427.2 


36 


J1837 


19121.227 


0.009 










J195900.31-252723.1 


46 


J1959 


20575.395* 


0.009 










J210318.76+021002.2 


16 


J2103 


19750.199 


0.009 










J210423.76+073140.4 


52 


J2104 


18065.091 


0.009 










J212454.61+203030.8 


21 


J2124 


19684.790 


0.009 


-0.41 


0.29 


-2.1 x 10" 5 


1 


J212808. 86+151622.0 


31 


J2128 


19426.310 


0.009 










J214510.25-494401.1 


18 


J2145 


19712.951 


0.009 










J220734.47+265528.6 


2 


J2207 


19978.751 


0.009 










J221058. 82+251 123.4 


45 


J2210 


18402.957 


0.009 










J224747.20-351849.3 


40 


J2247 


18853.775 


0.009 


+0.77 


0.64 


+4.1 x 10~ 5 


1 


J232607.07-294130.7 


8 


J2326 


19882.132 


0.004 










J234401. 81-212229.1 


44 


J2344 


18461.591 


0.007 


-0.313 


0.019 


-1.70 x 10- 5 


16 


J235333.60+455245.8 


7 


J2353 


19936.123 


0.009 











Notes. Period changes included only where significance is at least lcr. (n) In Norton et al. (201 1). <6) Period found here is longer than in Norton et al. 
d201lh . (c) GSC2314-0530. <rf) Data contaminated by nearby star's light curve. w ASAS Jl 11932-3950.8. (/) ROTSEI J170240. 1 1 + 15 1 122.7. 
<*> Insufficient data to calculate both O-C series. (,!) V1067 Her. (,) VI 104 Her. 



More notable is that three objects show a period change 
significant at more than 3 cr, where we would expect no sig- 
nificant results by chance alone in a set of 51 objects. These 
three stars underwent significant decreases in period during 
the time they were observed by SuperWASP: objects J 1743 



(-0.055 + 0.003 s yr 1 , significant at 21 cr), J 1331 (-0.075 + 
0.013 s yr" 1 , significant at 5 cr) and J2344 (-0.31 + 0.02 s yr 1 , 
significant at 16 cr). Figs. [T]-[6] show the folded light curves and 
O-C diagrams for these three objects, while Fig. Q illustrates 
a non-significant result (object J0509) for comparison. Given 
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Fig.3. Light curve for object J1331, folded at period of F i g . 5 . Light curve f or object J2344, folded at period of 
18836.380 s, with mean values for binned data overplotted. 18461.591 s, with mean values for binned data overplotted. 
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Fig. 4. O-C diagram for object J 1331, with best fit quadratic 
curves overplotted. Black stars and solid line indicate primary 
eclipse series; grey open diamonds and dashed line indicate sec- 
ondary eclipse series. Period decrease significant at 5 cr is indi- 
cated. 



Fig. 6. O-C diagram for object J2344, with best fit quadratic 
curves overplotted. Black stars and solid line indicate primary 
eclipse series; grey open diamonds and dashed line indicate sec- 
ondary eclipse series. Period decrease significant at 16 <x is indi- 
cated. 



the method of period determination, the choice of a precise zero 
point for the ephemerides of these three objects is somewhat ar- 
bitrary. However, using the zero points assigned by the program, 
and the optimum periods for the whole observation set, we have 
approximate eclipse times: 

J1743 HJD 2454997.688926 + 0.258108£ - 4.5 x 1CT I0 £ 2 
J1331 HJD 2454976.358606 + 0.218014£ - 5.2 x 10~ 10 E 2 
J2344 HJD 2454417.284057 + 0.213676£ - 2.1 x 10~ 9 £ 2 . 



4. Discussion 

We may note that peri od changes have been observed in many 
E Bs dKim et alJ 120031 cl aim that about 46% of the 1 140 EBs 
in iKreiner et alJ s Q2001) collection of O-C diagrams show at 
least some evidence of period change). Moreover, in short- 
period W UMa-type systems, period increases and decreases 
have been found to be of similar frequency, with the majority 
having \P\ < 0.02 s yr" 1 and none exceeding about +0.4 s yr~' 
(Kubiak et al. 2006), which would make our 3 main examples of 
period change notably rapid, but not exceptionally so. Therefore, 
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Fig. 7. O-C diagram for object J0509, with best fit quadratic 
curves overplotted. Black stars and solid line indicate primary 
eclipse series; grey open diamonds and dashed line indicate sec- 
ondary eclipse series. No significant period change is indicated. 



we do not suggest that our evidence for period change is in itself 
especially surprising, nor would we argue for period decreases 
being generally dominant in this stellar population, on the basis 
of such a small sample. Rather, the particular interest of (rela- 
tively rapid) period decrease in the three objects J 1743, J 1331 
and J2344 is that they are already near or below the short-period 
limit for W UMa objects of 0.22 d, and have timescales P/P be- 
tween 4 x 10 5 and 6 x 10 4 years. The underlying processes which 
may be operating in these objects therefore merit further consid- 
eration. 

To investigate possible causes of the period decreases ap- 
parently observed for these th ree objects, the eclip sing binary 
modelling s oftware PHOEBE dPrsa & Zwitterll2005l) . built upon 
the code of IWilson & Devinnevl d 197 lb . was used to estimate 
component masses, radii and other system parameters. In the 
absence of radial velocity measurements, values for semimajor 
axis a could not be determined, so best-fit parameter combina- 
tions were found for a range of physically-plausible values of 
a. Several external constraints were applied: from general light 
curve morphology, such as continuous light variation, the sys- 
tems were taken to be W UMa-type contact binaries (overcon- 
tact in PHOEBE's terms) i.e. with filling factors in (0,1]. The 
values for Mi were required to be approximately consistent with 
system temperatures as indicated by spectral types, themselves 
estimated from V — J and V - K colours using 2MASS data 
(J 1743 was taken to be approximately type K2, with tempera- 
ture ~ 4600 K; J 1331 type G9, r eff ~ 5000 K; J2344 type K9, 
T e ff ~ 3700 K). Mass ratios q were take n to be in the range 
0.08-0.8, as typical for W UMa systems (lHilditchll200"l . sub- 
ject to the hydrogen-burning requirement for M 2 to be at least 
O.O8M0. Maximum and minimum values of a could then be de- 
rived from the possible radii of the (contact) component stars. 
The resulting best-fit masses and radii (found to 2 s.f. by mini- 
mizing^ 2 values) are given in Table [2] Object J 1331 could not 
be as well-fitted as the other two objects, by any combination 
of a, q, inclination i and Kopal potential Q, because its maxima 
differ significantly in height, probably due to spots on the stellar 
surfaces. However, spots were not modelled here, in the absence 
of direct evidence. 



Period decrease might be caused by magnetic braking, re- 
moving angular mome ntum from the sy stem. The plausibi lity o f 
this was assessed using lBradstreet & Gurnard 's equation (Q 994): 



P a -1.1 x 10"V (1 +9) ( M i + M 2Y 



h 2 (M 1 R 4 l -M 2 R 4 2 )P'^(3) 



Although this appli es properly only to detached systems, 
Bradstreet & Guinan suggest that magnetic field strengths would 
be weaker in contact binaries. So we may at least use their for- 
mula to estimate an upper limit for the effect of magnetic brak- 
ing on a contact system. Therefore, taking t he gyration constant 
k 2 as 0. 1 (typical for main sequence stars, iBradstreet & Guinaru 
1994), the value of the RHS was evaluated for the various sys- 
tem parameter combinations and compared with dP/df for the 
three stars. The results, given in Table|2] are between about 20% 
and 1% of the observed quantities. Other estimates of the ef- 
fect of magnetic brakin g we re made from equati ons for J from 
Rappapor t et al.l (1 19831) and iHurlev et al.1 d2002l) . as calibrated 
bv lDavis et al.l (2008) by using the angular momentum loss rate 
at the upper edge of the cataclysmic variable period gap (the 
binary systems were approximated as single, fully-convective 
stars in order to use these equations). This produced significantly 
smaller estimates of expected period change: between 3 and 4 or- 
ders of magnitude too small. This suggests that magnetic braking 
is not the main cause of period decrease in any of the three sys- 
tems, though the difference of 1 to 2 orders of magnitude in the 
estimates produced by the different equations also implies sig- 
nificant limitations in our ability to quantify magnetic braking. 

Another possibility is angular momentum loss due to grav- 
itational wa ve r adiation (GWR) . Combining expressions from 
iKolbl d2010h and Hilditch] d200lh : 



g = 3^ = -1.27xl0-y MlM2 



(4) 







the expected P/P was calculated for each parameter combina- 
tion (Table However, these quantities are about 6 orders of 
magnitude smaller than the observed P/P values. GWR can be 
responsible for only a tiny fraction of the period decreases ob- 
served here. 

This leaves mass transfer from Mi to M 2 and/or mass (and 
hence angular momentum) loss from the system, from other 
mechani sms, as plausible causes. Using equations as given by 
iHilditchl : 



P 
P 



3Mi(M] -M 2 ) 



P 

~=3M, 



M\M 2 
(Mi + M 2 ) d 2 



M 2 



MiM 2 a 2 Mi(Mi+M 2 ) 



(5) 



(6) 



(where d is the distance from the binary centre of mass to 
Lagrange point L 2 ) for conservative and non-conservative mass 
loss respectively, we can calculate the necessary values of Mi 
to explain the observed P/P values (Table 0. These are of the 
order of 10~ 6 -10~ 7 Mq yr" 1 , which w e may note is similar to the 
mass loss rate calculated by Hilditch for contact binary S V Cen, 
to explain its precipitous observed rate of period decrease. He 
argues that mass may be ejected during a rapid phase of Roche- 
lobe overflow (or more generally, in a contact system) through 
L 2 , contributing to period decrease. Possible mechanisms for un- 
stable mass transfer/loss leading to rapid period d ecrease incon- 
tact or near-contact binaries are also discussed in lRasiol(IT99l . 
Tvle nda et al.l (1201 ll) and lJiang et alJ (1201 ll) . This explanation is 
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Table 2. Model parameters and theoretical period/mass changes for three objects. 
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left as the most plausible for the observed period decreases in the 
objects J 1743, J 1331 and J2344, with further small contributions 
from magnetic braking and GWR. 

What future, then, might we envisage for these three sys- 
tems? Since their primary masses and mass ratios cannot cur- 
rently be deter mined wit h any precision, i t is po ssible that they 
are subject to Jiang et alJ s low m ass limit (1201 ll) . le ading t o un- 
stable mass transfer, or to lRasiol 's tidal instabilities (1995), and 
will undergo merger on a relatively short timescale. iRasiol pro- 
poses the orbital decay time ?d ~ 10 3 — 10 4 years for an unstable 
W UMa system. A nother indicati on of the possi ble timescale for 
merger is given by Tvle nda et alj s study (1201 ll) of the decaying 
period of the contact binary progenitor of the V1309Sco out- 
burst. Using their exponential model for period decay, one can 
determine that their object would have been decreasing in period 
at a rate of 0.3 s yr" 1 about 130 years before the observed out- 
burst (and presumed merger). Since object J2344 appears to be 
currently undergoing period decrease at approximately this rate, 
and is already below the short-period limit, it is perhaps not in- 
conceivable that it might merge on a timescale of centuries or 
even decades from now (assuming, of course, that its period de- 
crease is maintained). Whatever their ultimate fate, these may 
be rare examples of objects caught in a brief transitional stage 
between stable states. 



5. Conclusions 

Our study of the periods of 53 W UMa candidate binary stars, 
observed with SuperWASP, confirmed that they are all very close 
to the short-period limit, and thus constitute a useful sample of 
sources for investigating the causes of this limit. In three of the 
objects, period decrease significant at 5 cr or more was found 
during their time of observation; the remaining objects' period 
changes were consistent with zero change at 3 cr. Modelling 
estimates for possible system parameters indicated that neither 
magnetic braking nor GWR are likely to be the primary cause 
of these period decreases. The rates of change observed in these 
three objects can best be explained by unstable mass transfer 
from primary to secondary components, and/or mass and angu- 
lar momentum loss from the systems, which could lead to merger 
on a relatively short timescale if the periods continue to decrease. 

These are potentially unusual transitional objects, capable of 
shedding light on the evolution of EBs and the causes of the 
short-period limit of W UMa binaries. As such they would re- 
pay further study and additional observations. In future we hope 
to obtain spectroscopic measurements of these systems, permit- 
ting more precise stellar modelling and parameter determina- 



tion. We also intend to run the programs developed here on other 
SuperWASP variable objects, to search for further interesting pe- 
riod variation. 
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